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Abstract 

The many-body system comprising a He nucleus, three electrons, and a positron has been studied 
using the exact diagonalization technique. The purpose has been to clarify to which extent the 
system can be considered as a distinguishable positronium (Ps) atom interacting with a He atom 
and, thereby, to pave the way to a practical atomistic modeling of Ps states and annihilation in 
matter. The maximum value of the distance between the positron and the nucleus is constrained 
and the Ps atom at different distances from the nucleus is identified from the electron and positron 
densities, as well as from the electron-positron distance and center-of-mass distributions. The 
polarization of the Ps atom increases as its distance from the nucleus decreases. A depletion of 
the He electron density, particularly large at low density values, has been observed. The ortho-Ps 
pick-off annihilation rate calculated as the overlap of the positron and the free He electron densities 
has to be corrected for the observed depletion, specially at large pores/voids. 
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I. INTRODUCTION 



The o-Ps atom is the bound state of an electron and a positron with total spin S — 1. In 
vacuum o-Ps has a relatively long lifetime of 142 ns because its annihilation via the fast two- 
gamma channel is prohibited by the conservation of the angular momemtum. Annihilation 
through the two-gamma channel with an electron of the matter having an opposite spin is 
possible when o-Ps interacts with matter. The resulting pick-off annihilation depends on the 
overlap of the o-Ps with the electron density of the matter and it can reduce the positron 
lifetime remarkably [1]. In metals and semiconductors the electron density is enhanced 
around the positron due to the attractive electron-positron interaction. This means that a 
Ps-like quasiparticle is created by the short-range screening of the positron by the electron 
cloud. In molecular matter and in some insulators such as Si02, the electron density is 
low in the interstitial regions and a Ps atom can be formed [2]. In this case the electron 
accompanying the positron hinders efficiently further screening of the positron by electrons 
in the matter. On the other hand, there is a finite overlap between the positron and the 
electrons of the material manifesting itself as the repulsive Pauli interaction and as the 
pick-off annihilation of o-Ps. 

The pick-off annihilation lifetime spectroscopy of orto-Positronium (o-Ps) has a rather 
unique role as a method capable to study open volumes in soft or porous materials in 
which Ps is formed. Size and density distributions of nanometre-scale voids in polymers [3] , 
porous Si02 [4, 5], and biostructures [6, 7] have been estimated by measuring the pick-off 
annihilation rate of o-Ps. However, in order to make a reliable quantitative analysis of 
the experimental results, one should be able to model o-Ps states and annihilation in these 
structures with predictive power. 

In contrast to an atom where the nuclear degrees of freedom can, to a good approximation, 
be treated classically, both the electron and the positron in the Ps atom are light quantum- 
mechanical particles and the non-adiabatic correlation effects have to be taken into account. 
For positrons in dense materials, such as metals and semiconductors, a good starting point is 
a separate calculation of the quantum- mechanical state of a single positron [8] . In practice, 
this calculation is based on the results of many-body theories for a delocalized positron in a 
homogeneous electron gas. The nature of the Ps atom would call for a quantum-mechanical 
many-body calculation of the interacting electron-positron system in the Coulomb field of 
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nuclei. In the case of materials systems, this is clearly beyond the present computer capacity. 
An attractive starting point for modeling Ps is to consider it as a distinguishable particle 
also when it is interacting closely with matter. This is in line with the view of a long-living 
o-Ps interacting with matter and annihilating through a pick-off process. 

The above-mentioned model of the Ps atom interacting as a single particle with matter 
finds support also from experimental and theoretical studies of Ps scattering off atoms and 
hght molecules. Recently, the scattering of Ps off noble gas atoms and hght molecules has 
been measured using a slow mono-energetic Ps beam [9]. The scattering properties of Ps 
have been seen to be dominated by the repulsive electron-electron Pauli interaction at short 
distances. At long distances the van der Waals attraction may play a role, especially in 
interactions with molecules of high polarizability. 

The positronium Hydride (HPs) bound state has been observed experimentally by 
Schader et al. [10] and its binding energy has been estimated to be around 1.1 eV (0.04 Ha). 
The bound states of Ps with several small atoms have been studied theoretically using Exact 
Diagonalization techniques. The most carefully studied system is HPs [11-13], which has 
been found to form a bound state with a binding energy of 0.039 Ha, in good agreement 
with the experimental value. LiPs [11, 14, 15] and NaPs [11, 15] systems have also been 
predicted to form bound states with binding energies of 0.0123 Ha and 0.0084 Ha, respec- 
tively. The calculations show that at large distances the positron overlaps with an extra 
electron forming a free Ps-like cluster, while, when the positron density overlaps with that 
of the electrons inside the atom the Ps-like character is lost. In the case of closed-shell noble 
gas atoms the strong electron-electron Pauli repulsion dominates the atom-Ps interaction 
and the formation of bound states is prevented. Thus far, only a small number of different 
atoms have been considered in this kind of accurate calculations. However, it is already 
clear that the atoms able to trap Ps have high electron affinities and open-shell electronic 
structures [16]. 

Mitroy and Ivanov have studied the scattering of o-Ps off rare-gas atoms and alkali 
metal ions. [17, 18] They determined the scattering lengths for wave vectors up to 0.5 Oq ^ 
corresponding to the Ps energy of ~1.7 eV. Besides rare-gas atoms, also molecules have 
small electron affinities. They typically have HOMO-LUMO gaps of up to 1 Ha and their 
electrons are well localized in molecular orbitals. In polymers and biomolecular materials, 
the main cohesive interaction between the molecules is the van der Waals interaction and 
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their molecular orbitals are not greatly affected. Individual molecules still show HOMO- 
LUMO gaps and the repulsive Pauli interaction between electrons in neighboring molecules 
plays a role. On the other hand, the van der Waals interaction between the molecules and 
Ps can also be strong, especially for molecules of high polarizability. 

A model widely used to study Ps in small pores and voids (< 2 nm) is that developed 
by Tao and Eldrup [19, 20]. The semi-empirical model relates the Ps lifetime and the 
pore/void radius. Ps is considered as a single particle and the pore/void is modeled using a 
spherical, cubic or elliptic infinite potential well. The electron density of the material enters 
the void and forms a layer of the thickness AR on its surface. The pick-off annihilation 
rate is calculated from the overlap of the positron density with the electron density of the 
material. The parameter of the model, AR, is not well known for all the materials and its 
chemical dependence cannot be described. In addition, the model assumes compact pores 
inside the material. The model has been extended to larger pores including the annihilation 
with excited particle-in-a-box states distributed according the Boltzmann distribution at a 
finite temperature. [21] 

Schmitz and Miiller-Plathe [22] introduced an atomistic model to describe Ps states 
and annihilation at voids in polymers. The potential landscape for Ps is based on the 
superposition of interactions between the Ps atom and individual atoms. These interactions, 
in turn, include the long-range van der Waals interaction based on the polarizabilities of the 
Ps atom and the atoms in the polymer and the short-range Pauli repulsion between the 
electron in the Ps atom and the electrons of the material. The latter interaction is described 
by fitting to experimental scattering cross sections. The Ps wave function corresponding 
to a given finite temperature is solved for by the path integral Monte-Carlo method and 
thereafter the pick-off annihilation rate is calculated as the overlap of the positron density 
with the electron density which is also obtained as a superposition over the atoms. 

Our final goal is the atomistic modeling of a Ps atom interacting with molecular materials 
so that the necessary approximations made and the model parameters chosen are based on 
ab-initio results. Studying theoretically the HePs system will help to understand several 
aspects in the properties of Ps also in molecular solids and liquids. This is because in 
molecules electrons are also in a closed shell configuration and a HOMO-LUMO gap sets the 
minimum energy necessary for an extra electron (of the Ps atom) to enter into the electron 
cloud of the molecule. 
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The similarities between molecules with HOMO-LUMO gaps and closed-shell rare gas 
atoms in mind, we have calculated the many-body wavefunction and the total energy of 
the unbound HePs system by the exact diagonalization technique using a explicitly corre- 
lated Gaussians basis optimized by a stochastic variational method (ECG-SVM). In order 
to study the interaction of a Ps and a He atom at finite distances, the mean distance of the 
positron from the nucleus has been constrained. As a result, a set of configurations has been 
considered with the nucleus-positron mean distance ranging between 1.87 ao and 91.90 ao. 
Electron and positron densities, interaction energies and annihilation rates have been cal- 
culated and analysed in order to clarify to which extent the picture of the distinguisable Ps 
particle can be applied. 



II. STOCHASTIC VARIATIONAL METHOD AND THE MODELING OF THE 
HePs SYSTEM 

We have described the A'"-particle system comprising the heavy nucleus (treated as a 
single particle), light electrons and the positron by the non-relativistic Hamiltonian 

where pj, mj, and q^, are the momenta, masses, and charges of the particles, respectively, 
is the distance between the i*'^ and -f^ particles, and Tcm is the center-of-mass (CM) kinetic 
energy. The wavefunction is written as a linear combination of properly antisymmetrized 
explicitly correlated Gaussian (EGG) [23] functions as 



* = E ^5M.(^, ^0 = 5^ Q ^ { exp (-1 YZ'^uIx ® ^SMs } , (2) 

i=i 1=1 

where A is an antisymmetrization operator and xsms is the spin eigenfunction with 
oxsMs ^ S{S + l)h'^XSMs and S^Xsms = MsfiXSMs- We have considered only the spher- 
ically symmetric total angular momentum L = states. The mixing coefficients Cj are 
obtained by diagonalizing the Hamiltonian matrix < '4^^sms\-^\'^^sms ^- '^^^ nonlinear co- 
efficients A^^j^ are randomly generated and a new ^4^^ is kept only if the update lowers the 
total energy of the system. Using the Jacobi coordinate set {xi,...,xjv-i}, with the reduced 
mass of the i^^ coordinate as — mj+i i^.i^l™'' , allows for a straightforward separation of the 



CM movement and therefore the total number of "particles" (degrees of freedom) treated is 
A^- 1. 

The stochastic variational method (SVM) is suitable for calculating many-body wave- 
functions, such as ^ in Eq. (2), for small systems. [23] The number of particles is limited 
by the present algorithms and computer capacity typically to less than seven particles. The 
ECG basis comprises in our calculations 600 functions. The matrix elements between the 
basis functions are computed analytically in an efficient manner. In order to keep the size of 
the ECG basis small, the basis set is optimized by choosing the most appropriate non-linear 
parameters A^,^ for the Gaussian functions using a stochastic search in the parameter space. 
The large number of parameters to optimize, up to tens of thousands, prevents the use of 
direct search methods. The best non-linear parameter for each pair is chosen after 25 trials 
within each optimization cycle and 20 optimization cycles are done for each basis function. 
The trial values are obtained choosing values between a certain minimum and maximum 
using a random number generator. 

SVM can be also used to study unbound systems setting boundaries to the possible values 
of the non-linear parameters A^^^, i. e. increasing the minimum value the particles can be pre- 
vented to separate to the infinity [17]. The configurations with the lowest energy are searched 
for within the available parameter space, i.e., the energy is minimized within a constrained 
parameter space. In order to model the interaction between He and the Ps atom at different 
distances it is enough to set a constraint only to the non-linear parameter ^4^^, corresponding 
to the nucleus-positron pair. Thus, the size of the parameter space available for optimizing 
the wavefunction is not seriously restricted. The nucleus-positron distance parameter Ihc-p 
corresponding to the maximum of the non-linear parameter, ^^e-p — (V^He-p)^' varied 
between 1 and 50 ao. The wavefunction and the total energy of the HePs system are then cal- 
culated for each constrained He-Ps distance. Due to practical reasons, the nucleus-electron 
distance parameter Ine-e also has to be limited. We have checked that using i^g-e — ^00 ao 
gives structures and energies well-enough converged for our discussion below. 

The kinetic energy ((T)) and the potential energy {{V)) of a system interacting through 
Coulomb potentials are related by the Virial theorem 2{T)—-{V). Any deviation from this 
relation is due to an inadequate basis function set. Optimizing the basis set the virial 
coefficient -|- 1) decreases in magnitude. 
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In order to study the many-body wavefunction of the system, the positron {pp{r) — 
(^|5(rp — rjv — and electron {pe{r) — X]^\(*|(^(^j,e — f^r — density distributions 
are plotted. Above, is the position vector of the nucleus and the index i runs over all the 
electrons A^e- The e-p pair correlation function is calculated as 

p,.,{f) = Y^imn - r> - (3) 

i=l 

The distribution of the electron-positron (e-p) CM is defined as 

Pe,''ir-) = J2{mr1'''-rW) (4) 

i=l 

where ff'^ = (ri,e + ^)/2 is the center of mass of the positron and the i^^ electron and i runs 
over all the Ne electrons in the HePs system. The calculated distributions do not have any 
angular dependency and the radial distributions are the solid-angle averaged distributions 
multiplied by the radius squared. 

The two electrons of the He atom are prepared in the S — state while the electron and 
the positron of the o-Ps are in the S = 1 state. The total spin angular momentum of the 
system is, then, S* = 1. For the pick-off annihilation rate r^°, the S* = component of the 
i*^electron-positron pair is selected through the PjIq ^Pi^ projector, i.e., 

= 47rrlcJ2mm - r,)Pn,\^) (5) 

i=l 

where Tq is the classical electron radius and c is the speed of light. The summation runs 
over all the electrons. 

III. RESULTS & DISCUSSION 

Before discussing the main topic of this work, the identification of the Ps atom within the 
HcPs system, we will shortly consider the energetics of the HcPs system. The interaction 
energy Ej of HePs is defined as the difference between the total energy of the interacting 
system and the sum of the energies Ehs and Eps of the isolated He and Ps atoms, respec- 
tively, i.e., Ei = EhcPs - Ehc - Ep^. £;ps=-0.24999999999 Ha and Ene = -2.903693749 Ha 
have been obtained from SVM calculations with a basis of 56 functions (virial coefficient = 
3.0x10"^^) and 600 functions (virial coefficient = 1.8x10""^°), respectively. The interaction 
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energy is shown in figure 1 as a function of the nucleus-positron mean distance (rp). Ej is 
positive and approaches zero for large (r^) values. For {vp) smaller than 2.5 ao, Ej is larger 
than Eps. Our calculations do not show the formation of the bound He-e^ system. He can 
bind a positron only when it is in the S = 3 state [24] and its energy is 0.4 Ha above the the 
highest energy we have considered for the HePs system. It should be noticed that, already 
at 5 ao, the interaction energy is clearly above the energy of a thermalized Ps. 
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Figure 1. Total interaction energy and its kinetic and potential energy contributions as a function 
of {rp) . The horizontal dotted line denotes the energy above which the Ps atom is unstable against 
dissociation. The inset shows the log-log plot of the virial coefficient versus (rp). 



The error of the interaction energy AE has been estimated as the total energy decrease 
after an optimization loop of the function basis. For the optimization loop the best non-linear 
parameter has been chosen for each pair after 25 trials. This procedure has been repeated 
20 times for each basis function. For the most confined systems, {rp)= 2 ao, AE = 0.01 Ha 
and the virial coefficient is ~ 0.1. The confinement method shrinks the parameter space 
available to optimize the basis. In addition, the many-body wavefunction has to include the 
strong correlations happening in confined systems. Accordingly, the convergence is slow and 
the virial coefficient is large. For (vp) = 48 ao, the error and the virial coefficient decrease 
to AE = 2x10"^ Ha and 7x10"^, respectively. 

Figure 1 also shows the decomposition of the interaction energy to the kinetic and poten- 
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tial energy contributions. They are defined similarly to the total interaction energy. F^xin 
is the main contribution to the interaction energy at distances shorter than 10 ao, mostly 
due to the confinement of Ps and its small mass. Epot is always positive and at very short 
distances (< 4 ao) it grows very fast. 

The calculated mean distances of electrons (rg), the positron (r^) and the electron- 
positron CM (r^'^p) from the He nucleus are shown in Fig. 2 as a function of the nucleus- 
positron distance parameter Ine-p corresponding to the confinement of the system. The 
mean distances grow monotonously because HePs is ultimately unbound. The fact that (rg) 
increases with (r^) is an indication of the formation of Ps. At Ine-p < 5 ao, (rg) saturates to 
a value of 2 ao because the electron of Ps does not penetrate into the He core. On the aver- 
age, it stays at ~ 4 ao from the nucleus (see figure 4). On the other hand, (r^) continues 
to decrease because the positron keeps being confined. 
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Figure 2. (Color online) Positron (red squares), electron (blue triangles) and electron-positron CM 
(black circles) mean distances from the He nucleus for different calculated configurations versus 
the nucleus-positron distance parameter Ine-p- The electron mean distance (blue line) obtained 
from the unpolarized Ps model is also given. 

Using (vp) and assuming that the Ps atom is not polarized, (rg) can be approximated [25] 
as the sum of two thirds of the mean electron distance in the He atom of 0.93 ao and one 
third of the mean positron distance (r^) (blue line in figure 2). The agreement with the (rg) 
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calculated from the many-body wavefunction is very good when (r^) is larger than 7-8 ao 
(or the nucleus-positron distance parameter Ine-p is larger than 4-5 ao). For smaller Ihc-p 
values the actual (ve) is larger than the estimated one which reflects the polarization of 
the Ps atom when approaching the He nucleus. The difference can actually be taken as a 
measure of the polarization (see figure 3). For (rp)=20 ao the polarization parameter is very 
small (~6xl0~^ ao). The polarization starts to increase when (r^) is below 10 ao and at 
4 ao its value is ~0.2 ao , i.e., 10% of the electron-positron mean distance for an isolated Ps. 
The maximum value of the polarization parameter is 0.8 ao when the interaction energy is 
0.44 Ha and Ps is actually unstable against dissociation. According to our calculations, the 
smallest stable configuration corresponds to {vp) ~ 2.9 ao and a polarization parameter of 
~0.45 ao. 




Nucleus-positron mean distance (ao) 



Figure 3. The polarization parameter for all the calculated HePs configurations as a function of the 
positron mean distance (vp) from the He nucleus. The inset shows the electron-positron contact 
density corresponding to the Ps atom component as a function of (r^) . The solid Une indicates the 
contact density of free Ps. 



The electron-positron contact density of the Ps atom (shown in the inset of figure 3) has 
been obtained by subtracting the contribution of the He electrons estimated as two times 
the electron-positron contact density of the S—0 electron-positron pair. At (rp)=20 ao, the 
contact density within the Ps atom reaches the value of 0.97, which is close to the exact 
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value of 1 for an isolated Ps. At (rp)=4 ao the contact density is still ~0.94 and at the 
smallest separation decreases down to 0.76. The decrease of the contact density reflects and 
increasing polarization of the Ps atom in the interaction with the He atom. 

Figure 4 shows the radial distributions of the positron and the electrons in HePs config- 
urations corresponding to the positron mean distance (rp) ranging from 3.4 ao to 13.7 ao. 
The electron distribution can be decomposed into two components. The larger corresponds 
closely to the undisturbed electron density of a He atom. The second maximum overlaps 
with the positron distribution especially well for large (rp) indicating the formation of a 
well-distinguishable Ps atom. At small separations, the overlap of the electron and positron 
distributions weakens, indicating the increase of the polarization of the Ps atom. 
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Figure 4. (Color online) Radial distributions of the positron (full lines) and electrons (dashed lines) 
for four different configurations of the HePs system with (rp) values of 3.1 ao (light green), 5.2 ao 
(red), 7.5 ao (dark green) and 12.1 ao (blue). The arrowed tags indicate the position and the value 
of the positron mean distances (rp) from the He nucleus for each configuration. The shading gives 
the electron density of a free He atom. All the distributions are normalized to the total number of 
particles. 

Figure 5 shows the radial electron-positron correlation distribution for the four HePs 
configurations and for the isolated Ps atom. The distribution is calculated from equation (3) 
and it is normalized to the total number of electrons. When (rp) = 3.4 ao, the electron- 
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positron pair correlation function shows a single peak but when (rp) is above 5 ao a second 
peak starts to appear, as a shoulder, at 2 ao- The shape and the magnitude of the second 
peak approach clearly to those for the isolated Ps atom and it corresponds to the electron- 
positron correlation inside the Ps atom. The second peak follows the positron distribution 
when it recedes from the He nucleus and it corresponds to the correlation between the 
positron and the electrons of the He atom. However, even at 8.5 ao, the interaction between 
He and Ps affects the electron-positron correlations and the appearance of the Ps atom is 
not so clear-cut as in the electron distribution in figure 4. 
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Figure 5. Electron-positron pair correlation function for the four HePs configurations of figure 4. 
The correlation distance for the isolated Ps atom (the dashed filled curve) is shown for comparison. 
All the electron-positron correlation curves are normalized to the number of electrons in the system 
in question. 



The distribution of the e-p CM, shown in figure 6, can also be decomposed into two 
components. The larger component is closer to the He nucleus and its integrated intensity 
is double of that for the smaller and more delocalized component. Thus, the larger and 
the smaller components correspond to the two electrons of the He atom and the electron of 
the Ps atom, respectively. The Ps component is very similar to the positron radial density, 
especially when the mean He nucleus-positron distance is larger than 5 ao. 

Figure 7 shows the positron pick-off annihilation rate P^" as a function of (r^), calculated 
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Figure 6. (Color online) Electron-positron CM distribution for the four HePs configurations. The 
electron (dashed) and the positron (filled areas) radial distributions are also shown for comparison. 
The CM distributions are normalized to the number of electrons. 

using equation 5. It is an important magnitude for understanding Ps lifetime experiments 
in molecular matter and it also reflects the interaction of the Ps and He atoms and its 
behaviour. The projector PgZ^Q in equation 5 ensures that the pick-off processes occur only 
with the He electron forming a spin-singlet state with the positron. F^" increases with the 
increasing overlap of the positron and He electron density. At (rp)=10 ao, F^" is about the 
same as the self-annihilation rate ~7.04x 10"'^ 1/ns of o-Ps (lifetime r ~ 142 ns). At 4 ag, 
TP" ~ 0.2 1/ns (r ~ 4.8 ns) and for (rp) ~ 2.5 a^ Fp° equals to 2 1/ns (r=500 ps), the 
spin-averaged annihilation rate of Ps. 

A non-selfconsistent positron pick-off annihilation rate is obtained from the overlap of 
the free He electron density n^J and the SVM positron density as 



The free He electron density is from a SVM calculation and the factor one half reflects that 
only one electron in He annihilates with the positron. F^° is widely used to calculate the 
pick-off annihilation rate of o-Ps [26]. With increasing nucleus-positron mean distance, it 
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Figure 7. (Color online) Positron annihilation rate obtained from the SVM results using equation 5 
(black circles) and from the overlap of the free He atom and the positron densities, using equation 6 
(solid line). The inset show the ratio of the values. The blue horizontal line gives the self- 
annihilation rate of o-Ps. 

agrees with the many-body result until (r'p)~5 ao (solid line in figure 7). When (r^) 
further increases, r*^° remains higher than T^". The ratio r^°/rP°, shown in the inset of 
figure 7, ranges between 1 and 0.6. The most confined systems show the largest ratios and 
it decreases steadily until (vp) ~ 30 ao. At the level of the o-Ps self-annihilation rate the 
overlap of the annihilating positron with the He electrons has decreased by 20% compared 
to the undistorted He electron density. It should be noted that at large separations the 
overlap and the annihilation rate are very small, which is also reflected in the scatter of the 
overlap ratios. However, small values down to o-Ps self-annihilation rate are expected to be 
important for large voids. 

In contrast to a positron in an electron gas, where the positron-electron contact density is 
enhanced [8] , for Ps interacting with He the many-body effects induce a depletion of the con- 
tact density. The electron in Ps screens the charge of the positron and the electron-electron 
Pauli repulsion further decreases the contact density at the positron position, specially at 
low electron density values. Thus, in strongly confined systems the contact density is rel- 
atively weakly affected, but in weak confinement it is strongly depleted. needs to be 
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corrected for this effect specially when calculating the annihilation rate in large voids or 
pores, where the contact density is small and the annihilation rate is largely overestimated. 

Figure 8 resumes the polarization parameter, He nucleus-positron mean distance and the 
interaction energy of the HePs system as a function of the positron lifetime (l/r=r^°+l/142.05). 
The gray area marks where the interaction energy of Ps is larger than the Ps binding energy. 
The figure gives a clear idea about the magnitude of these quantities in different conditions. 



CD 

>■ 

a> 


c 

c 
o 

o 



o 
Q. 



o 



0.5 
0.4 

0.3 HO 
0.2 
0.1 



<> 



0- 



8 
6 
4 

2H: 



0.8 
0.6- 
0.4- 
0.2 

0- 



Ps formation 



Oo 



C> , O O .O 



□ □ 



Ps formation 



-3 



O 
O 



Ps formation 



O O o o O o 



T I I I I I I I I I r 

10 20 30 40 50 

Lifetime (ns) 



OOO 



60 



Figure 8. (Color online) The polarization (black circles), the He nucleus-positron distance (blue 
squares) and the interaction energy (green diamonds) of the HePs system as a function of the 
positron lifetime. The grey area indicates the value range where the interaction energy of Ps atom 
is larger than its binding energy. 



IV. CONCLUSIONS 

We have been interested in to what extent a separate Ps atom can be distinguished in a 
close interaction with a He atom, i.e., when the volume of the HePs system is constrained. 
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According to our many-body calculations the behavior of several measures, i.e., those of 
electron and positron mean distances from the He nucleus, electron and positron densities, 
as well as electron-positron pair correlation function and center-of-mass distribution show 
the formation of a Ps atom that is weakly polarized until the mean He nucleus-positron 
mean distance (r^) is shorter than ~5 ao- At shorter distances the Ps atom also forms but 
it is strongly polarized. The most clear-cut evidence is given, simply, by the electron and 
positron densities shown figure 4. There, the electron density can be decomposed into the He 
and Ps atom components and the positron density overlaps well with the electron density. 
Even at (r^) fa 5 ao the positron-electron contact density is still high, ~93 % of that in the 
free Ps atom, and the polarization parameter indicating the imbalance of the electron and 
positron densities in the Ps atom is small, 0.2 ao. 

The above notions for the HePs system, generalized to molecules, supports the possibihty 
to use a practical single-particle description for o-Ps in molecular materials. Many-body 
calculations could be used to define pairwise atom-Ps potentials which could form the basis 
for the calculation of the Ps energy landscape in the material and obtain, subsequently, its 
distribution and pick-off annihilation lifetime. In that respect, the above-mentioned findings 
about the relative depletion of the contact density, will simplify the reliable estimation of 
the pick-off annihilation rate. 

ACKNOWLEDGMENTS 

This work was supported by the Academy of Finland through the postdoctoral researcher 
fellowship and the centre of excellence programs. Thanks are due to K. Varga for proving us 
his the ECG-SVM code. We also thank 1. Makkoncn, J. Mitroy and T. Rantala for helpful 
discussions. The comments and suggestions from I. Makkonen have been of great help while 
preparing the manuscript. 



[1] O. E. Mogensen, in Positron Annihilation in Chemistry, Vol. 58, edited by H. K. V. Lotsch 
(Springer- Verlag, 1995) Chap. Principles of Positron Annihilation in Molecular Solids, p. 193. 
[2] S. V. Stepanov and V. M. Byakov, Principles and application of positron & positronium chem- 

16 



istry, edited by Y. C. Jean, P. E. Mallon, and D. M. Schrader, 5 (World Scientific, 2003) Chap. 

Physical and radiation Chemistry of the positron and positronium, p. 117. 
[3] A. Uedono, R. Suzuki, T. Ohdaira, T. Uozumi, M. Ban, M. Kyoto, S. Tanigawa, and 

T. Mikado, J. Polym. Sci. Part B 36, 2597 (1998). 
[4] Y. Nagai, Y. Nagashima, and T. Hyodo, Phys. Rev. B 60, 7677 (1999). 
[5] L. Liszkay, C. Corbel, P. Perez, P. Desgardin, M.-F. Barthe, T. Ohdaira, R. Suzuki, P. Crivelh, 

U. Gendotti, A. Rubbia, M. Etienne, and A. Walcarius, Appl. Phys. Let. 92, 063114 (2008). 
[6] P. Sane, E. Salonen, E. Falck, J. Repakova, F. Tuomisto, J. Holopainen, and I. Vattulainen, 

J. Phys. Chem. B letters 113, 1810 (2009). 
[7] A. W. Dong, C. Pascual-Izarra, S. J. Pas, A. J. Hill, B. J. Boyd, and C. J. Drummond, J. 

Phys. Chem. B 113, 84 (2009). 
[8] M. J. Puska and R. M. Nieminen, Rev. Mod. Phys. 66, 841 (1994). 

[9] S. J. Brawley, S. Armitage, J. Beale, D. Leslie, A. Williams, and G. Laricchia, Science 330, 
789 (2010). 

[10] D. M. Schrader, F. M. Jacobsen, N.-P. Frandsen, and U. Mikkelsen, Phys. Rev. Lett. 69, 57 
(1992). 

[11] G. G. Ryzhikh, J. Mitroy, and K. Varga, J. Phys. B: At. Mol. Opt. Phys. 31, 3965 (1998). 

[12] J. Mitroy, Phys. Rev. A 73, 054502 (2006). 

[13] S. Bubin and K. Varga, Phys. Rev. A 84, 012509 (2011). 

[14] J. Mitroy, J. At. Mol. Sci. 1, 275 (2010). 

[15] J. Mitroy and G. G. Ryzhikh, J. Phys. B: At. Mol. Opt. Phys. 34, 2001 (2001). 

[16] J. Mitroy, M. W. J. Bromley, and G. G. Ryzhikh, J. Phys. B: At. Mol. Opt. Phys. 35, R81 

(2002). 

[17] J. Mitroy and I. A. Ivanov, Phys. Rev. A 65, 012509 (2001). 

[18] J. Mitroy, J. Y. Zhang, and K. Varga, Phys. Rev. Lett. 101, 123201 (2008). 

[19] S. J. Tao, J. Chem. Phys. 56, 5499 (1972). 

[20] M. Eldrup, D. Lightbody, and J. N. Sherwood, Chem. Phys. 63, 51 (1981). 
[21] T. L. Dull, W. E. Frieze, D. W. Gidley, J. N. Sun, and A. F. Yee, J. Phys. Chem. B 105, 
4657 (2001). 

[22] H. Schmitz and F. Miiller-Plathe, J. Chem. Phys. 112, 1040 (2000). 
[23] K. Varga and Y. Suzuki, Phys. Rev. C 52, 2885 (1995). 

17 



[24] J. Mitroy, Phys. Rev. A 72, 032503 (2005). 

[25] A. Zubiaga, F. Tuomisto, and M. Puska, "Study of unbound HePs using exact diagonalization 

technique," (2012), to be published in the Material Science Forum. 
[26] W. Brandt, S. Berko, and W. W. Walker, Phys. Rev. 120, 1289 (1960). 



18 



